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I Abstract 

o . 

, We solve the problem of formulating Brownian motion in a relativis- 

ts) ■ tically covariant framework in 3 + 1 dimensions. We obtain covariant 

Fokkcr-Planck equations with (for the isotropic case) a differential oper- 
C*") , ator of invariant d'Alembert form. Treating the spacelike and timelike 

fluctuations separately in order to maintain the covariance property, we 
| show that it is essential to take into account the analytic continuation 

Q j. of "unphysical" fluctuations. We discuss the notion of locality in this 

I 1 framework and possible implications for entangled states. 

-I— > 
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1 Introduction 

The ideas of stochastic processes originated in the second half of the 19th century 
in thermodynamics, through the manifestation of the kinetic theory of gases. 
In 1905 A. Einstein £Q in his paper on Brownian motion provided a decisive 
breakthrough in the understanding of the phenomena. Moreover, it was a proof 
convincing physicists of the reality of atoms and molecules, the motivation for 
Einstein's work. It is interesting to note that Einstein predicted the so called 
Brownian motion of suspended microscopic particles not knowing that R. Brown 
first discovered it in 1827 0- The resemblance of the Schrodinger equation to the 
diffusion equation had lead physicists (including Einstein and Schrodinger) to 
attempt to connect quantum mechanics with an underlying stochastic process, 
the Brownian process. 

Nelson in 1966, constructed the Schrodinger equation from an analysis 
of Brownian motion by identifying the forward and backward average velocities 
of a Brownian particle with the real and imaginary parts of a wave function. He 
pointed out that the basic process involved is defined non-relativistically, and 
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can be used if relativistic effects can be "safely" neglected. The development 
of a relativistically covariant formulation of Brownian motion could therefore 
provide some insight into the structure of a relativistic quantum theory. 

Nelson pointed out that the formulation of his stochastic mechanics in the 
context of general relativity is an important open question 0] . The Riemannian 
metric spaces one can achieve, in principle, which arise due to nontrivial cor- 
relations between fluctuations in spatial directions, could, in the framework of 
a covariant theory of Brownian motion, lead to spacetime pseudo-Riemannian 
metrics in the structure of diffusion and Schrodinger equations. Morato and 
Viola have recently constructed a relativistic quantum equation for a free 
scalar field. They assumed the existence of a 3D (spatial) diffusion in a co- 
moving frame, a non-inertial frame in which the average velocity field of the 
Brownian particle (current velocity) is zero. In this frame the location of the 
Brownian particle in space experiences Brownian fluctuations parametrized by 
the proper time of the comoving observer. They interpreted possible negative 
0-component current velocities with what they called 'rare events', which are 
time reversed Brownian processes (a peculiarity arising in the relativistic treat- 
ment). The equation they achieved this way is approximately the Klein-Gordon 
equation. It is important to note that in the inertial frame they do not obtain a 
normal diffusion. This is due to the fact that their process is stochastic only in 
three degrees of freedom and therefore is not covariant. In this paper we shall 
study a manifestly covariant form of Brownian motion. 

Nelson himself 4 argued that Markov processes may lead to inconsistencies 
(which he demonstrated in entangled systems such as E.P.R |H])unless a non- 
local theory is allowed, therefore according to his own preference for a local 
theory he suggested to consider non-Markovian theories. We shall treat this 
problem in our concluding section. 

In a previous work |7| we introduced a new approach to the formulation of 
relativistic Brownian motion in 1+1 dimensions. The process we formulate is 
a straightforward generalization of the standard one dimensional diffusion to 
1+1 dimensions (where the actual random process is thought of as a 'diffusion' 
in the time direction as well as in space), in an inertial frame. The equation 
achieved is an exact Klein-Gordon equation. It is a relativistic generalization of 
Nelson's Brownian process, the Newtonian diffusion. In this paper we review the 
relativistic Brownian process in 1+1 dimensions [x| where the inclusion of both 
spacelike and timelike motion for the Brownian particle (event) is considered; 
if the timelike motion is considered as "physical" the "unphysical" spacelike 
motion is represented (through analytic continuation) by imaginary quantities. 
We extend the treatment to 3+1 dimensions using appropriate weights for the 
imaginary representations. The extension of the process to a general covariant 
form will be carried out in a succeeding paper. The complete formalism then 
can be used to construct relativistic general covariant diffusion and Schrodinger 
equations with pseudo-Riemannian metrics which follow from the existence of 
nontrivial correlations between the coordinate random variables. 

Finally, we discuss the possible implications of the process we consider (i.e. 
a relativistic stochastic process with Markov property which preserves macro- 
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scopic Lorentz covariance) on the entangled system, where we claim that though 
fluctuations which exceed the velocity of light occur the macroscopic behavior 
dictated by the resulting Fokker-Planck equation is local. 

2 The Problem of Assigning a Dynamical Evolu- 
tion Parameter for the Relativistic Brownian 
Process and the Stueckelberg Formalism 



Brownian motion, thought of as a series of "jumps" of a particle along its path, 
necessarily involves an ordered sequence. In the nonrelativistic theory, this 
ordering is naturally provided by the Newtonian time parameter. In a relativistic 
framework, the Einstein time t does not provide a suitable parameter. If we 
contemplate jumps in spacetime, to accommodate a covariant formulation, a 
possible spacelike interval between two jumps may appear in two orderings in 
different Lorentz frames. The introduction of proper time as a parameter for 
the RBP (Relativistic Brownian Process) is not adequate since in this case the 
second order correlations in the simplest case (i.e. for an isotropic homogeneous 
process with a diffusion constant a 2 ) have the form: 

E^Ax^Ax^) = 2a 2 As (1) 

for each [i] however, summing over /x, 

E(AxpAx?) = As 2 oc As, (2) 

where the first equality is by the definition of proper time and the second equal- 
ity is due to the Brownian property expressed in Eq. (JTJ; there is an obvious 
contradiction. We therefore adopt the invariant parameter r as the dynamical 
variable for the Brownian process, first suggested in 1941 by E.C.G Stueckel- 
berg [H] • The introduction of the notion of an invariant time r permitted the 
discussion of world lines not monotonic in the ordinary (Einstein) time, t. On 
such world lines two points may have the same t coordinate and represent the 
occurrence of two particles existing at the same time. The time coordinate, t, 
may increase or decrease as r evolves (just as a particle in Newtonian mechanics 
may move in the positive or negative spatial directions with time) representing 
particles and antiparticles respectively. In this way Stueckelberg was able to 
describe pair-creation and annihilation on a classical level and write a relativis- 
tic Schrodinger equation (with four coordinates of spacetime and the invariant 
time as a parameter). 

In 1973 Horwitz and Piron (jj] developed this concept, suggesting that all 
physical systems evolve through one universal invariant time, not affected by 
interaction or dynamical coordinate transformations, just as in the Newtonian 
theory. The theory formally resembles Newton's theory with the most significant 
difference that Euclidian space is replaced by Minkowski spacetime. 
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The point particle, in Newtonian mechanics, represented by its spatial co- 
ordinates at each given time, is replaced by a new concept, the event; it is an 
object, evolving in spacetime, represented by four coordinates at each given 
(universal) time r. The collection of the event's (four) positions in spacetime at 
all (universal) times forms the worldline. The time-coordinate of an event in a 
given inertial frame is the time in which, according to the clocks of this frame 
the event is detected, though it could have occurred, in principle, at a different 
(universal) time, r. In particular let us consider a parabolic worldline (figure 
1), describing pair annihilation, at a given time toB, such that the line t — t± 
intersects the worldline twice. At the time t\ two events will be detected on the 
frame's clocks, one corresponding to a particle and the other to an antiparticle, 
each occurring at different (universal) time t. The interpretation of the event 
going backwards in time as the antiparticle was given first by Stueckelberg and 
was later used by Feynman; it is now an accepted concept. However, in the 
parabolic worldline there are segments in which the event goes faster than light 
(either forward or backward in t). Do these segments have any physical rep- 
resentation? Should they be included in a physical theory? As we shall see in 
this paper the answer for the second question is positive in the formulation of 
our relativistic Brownian process, and the occurrence of this state of motion 
is dictated by the demand of achieving Lorentz invariance (more explicitly the 
d'Alembert) operator in the relativistic diffusion equation. 



figure 1 : Particles and Events. A particle/antiparticle corresponds to a 
world line segment formed of a trajectory of an event whose time coordinate is 
monotonically increasing/decreasing with r. At t\ a particle and an anti-particle 
may be detected, both generated by an event B. The particle/antiparticle corre- 
sponds to the (left/right) branch of the parabola (within the forward/backward 
light cones) where the time coordinate is increasing/decreasing with r. The 
segment of motion outside of the lightcones are tachyonic (spacelike) and are 
required by continuity of the worldline. These two points occur at different r 
but at the same t. At times t > toB (the turning point), the event B does 
not exist therefore the worldline generated by event B corresponds to a pair 
annihilation occurring at t = toB- An event, say, C may generate a worldline 
corresponding to a pair creation occurring at t = toe and another event A may 
generate a monotonic worldline therefore corresponding to a single particle. 

3 The Negative Correlation Problem 

A second fundamental difficulty in formulating a covariant theory of Brownian 
motion lies in the form of the correlation function of the random variables of 
spacetime. The correlation function for the isotropic (non-relativistic) Wiener 
[TU] process, is given by, 

< Aa^Aa^ >= a l ' J 'dt 
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a 1 ! = 2a 2 6 ij 

i,j = 1,2,3 (3) 

where the 5^ generates the Euclidian structure for the manifold on which the 
Brownian evolution take place The straightforward covariant generalization 
to the relativistic case is, 

<dw,(r)dw u (r')>= 2a2 ^ dT itr' W 

where rf v = diag(— 1, 1, 1, 1) is the Minkowski metric, and therefore 
< dwo(T)dwo{T) > < 0, which is impossible. Let us consider, however, a pro- 
cess which is physically restricted to only to spacelike or timelike jumps. One 
may argue that Brownian motion in spacetime should be a generalization of the 
non-relativistic problem, constructed by observing the non-relativistic process 
from a moving frame according to the transformation laws of special relativity. 
Hence the process taking place in space in the non-relativistic theory would be 
replaced by a spacetime process in which the Brownian jumps are spacelike. 
The pure time (negative) self-correlation does therefore not occur. In order to 
meet this requirement, we shall use a coordinatization in terms of generalized 
polar coordinates which assures that all jumps are spacelike. Consider for ex- 

2 

ample a relativistic Brownian probability density of the form e _ ^37, where [i 
is the invariant spacelike interval of the jump. This is a straightforward gen- 
eralization of the standard Brownian process in 3D, which is generated by a 

probability density of the form e^r , where r is the rotation invariant (i.e. the 
vector length)and a is proportional to the diffusion constant. We shall refer to 
this function as the relativistic Gaussian. As we shall see, a Brownian motion 
based on purely spacelike jumps does not, however, yield the correct form for an 
invariant diffusion process. We must therefore consider the possibility as well 
that, in the framework of relativistic dynamics, there are timelike jumps. The 

CT 2 

corresponding distribution would be expected to be of the form e~^7, where 
a is the invariant interval for the timelike jumps, and b is some constant. By 
suitably weighting the occurrence of the spacelike process (which we take for our 
main discussion to be "physical" , since its nonrelativistic limit coincides with the 
usual Brownian motion) and an analytic continuation of the timelike process, 
we show that one indeed obtains a Lorentz invariant Fokker-Planck equation 
in which the d'Alembert operator appears in place of the Laplace operator of 
the 3D Fokker-Planck equation. One may, alternatively, consider the timelike 
process as "physical" (as might emerge from a microscopic model with scatter- 
ing) and analytically continue the spacelike ( "unphysical" ) process to achieve a 
d'Alembert operator with opposite sign. 
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4 Brownian motion in 1+1 dimensions 



We consider a Brownian path in 1 + 1 dimensions generated by a stochastic 
differential (analogue to the Langevin equation ^2] an d Smoluchowsky process 
of the form: 

dx»(r) = ^{x{r))dT + dw^{r) 1 (5) 

where dw is a random process, which is a relativistic generalization of the Wiener 
process, whose properties will be defined later, and /3 M is a deterministic field 
(the drift). 

We start by considering the second order term in the series expansion of a 
function of position of the particle on the world line, /(x M (r) + Ax^), involving 
the operator 

= Ax fi Ax v a ° . (6) 
dx^ dx v 

We have remarked that one of the difficulties in describing Brownian motion in 
spacetime is the possible occurrence of a negative value for the second moment 
of some component of the Lorentz four vector random variable. If the Brownian 
jump is timelike, or spacelike, however, the components of the four vector are 
not independent, but must satisfy the timelike or spacelike constraint. Such 
constraints can be realized by using parameterizations for the jumps in which 
they are restricted geometrically to be timelike or spacelike. We now separate 
the random jumps into space-like jumps and time-like jumps accordingly, i.e., 
for the spacelike jumps, 

Aw 1 = i/icosha , Aw° = ^sinha (7) 

and for the timelike jumps, 

Aw 1 = a sinliQ , Aw° = iercosha (8) 

Here we assume that the two sectors have the same distribution on the hyper- 
bolic variable. We furthermore assume that /z,cr are generated by a relativistic 
Gaussian distribution, working in a Lorentz frame where the a distribution is 
assumed to be independent of /x, a and is uniformly distributed on the restricted 
interval [-L, L] (see discussion below) where L is arbitrary large. Therefore, in 
this frame < Aw'" > is (this is true in all frames; see discussion in Section 5) 
and we pick a normalization such that (for any component) < Aw n >oc At? so 
to first order in At the contribution to < O > comes only from < Aw^Aio" >. 

For a particle experiencing space-like jumps only, the operator O takes the 
following form: 



Q2 q1 q2 

''spacelike = V? [cosh 2 a ^-^ + 2 sinh a cosh a -Q-Q- t + sinn2 a Q^^ 
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If the particle undergoes time-like jumps only the operator O takes the form: 



d 2 d 2 d 2 

Oumeiike = cr 2 [sinh 2 a ^-^ + 2 sinh a cosh a g^gi + cosh 2 a-^-] (10) 

Since fi,cr and a are random processes, the average value of the operator O 
is the sum of the two averages of Eq. (JHJ and Eq. HlOjl . A difference between 
these two averages, leading to the d'Alembertian operator can only be obtained 
by considering the analytic continuation of the timelike process to the spacelike 
domain, choosing /i 2 = — a 2 . 

This procedure is analogous to the effect, well-known in relativistic quan- 
tum scattering theory, of a physical process in the crossed (t) channel on the 
observed process in the direct (s) channel. For example, in the LSZ formulation 
of relativistic scattering in quantum field theory [e.g. HH], a creation operator 
in the "in" state may be moved to the left in the vacuum expectation value ex- 
pression for the S'-matrix, and an annihilation operator for the "out" state may 
be moved to the right. The resulting amplitude, identical to the original one 
in value, represents a process that is unphysical; its total "energy" (the sum of 
four-momenta squared) now has the wrong sign. Assuming that the S'-matrix is 
an analytic function, one may then analytically continue the energy-momentum 
variables to obtain the correct sign for the physical process in the new channel. 
Although we are dealing with an apparently classical process, as Nelson has 
shown, the Brownian motion problem gives rise to a Schrodingcr equation, and 
therefore contains properties of the differential equations of the quantum theory. 
We thus see the remarkable fact that one must take into account the physical 
effect of the analytic continuation of processes occurring in a non-physical, in 
this case timelike, domain, on the total observed behavior of the system. 

In the timelike case, the velocity of the particle Aw 1 / AnP < 1. We shall 
here use the dynamical association of coordinate increments with energy and 
momentum 



so that 



E = M— P = M—, (11) 

a 2 = {^)\E 2 - P % (12) 

where M is a parameter of dimension mass associated with the Brownian par- 
ticle. It then follows that E 2 — p 2 = (-^f) 2 ^ 2 > 0. For the spacelike case, 
where p/E > 1, we may consider the transformation to an imaginary represen- 
tation E -> iE' and p -> ip', for E',p' real 1 , but E 2 - p 2 -»• p' 2 - E' 2 > 0. 
In this case, we take the analytic continuation such that the magnitude of a 2 
remains unchanged, but can be called — /i 2 , so that E' 2 — p' 2 = ^ 2 with n imag- 
inary. The spacelike contributions are therefore obtained in this mapping by 
E,p — > iE,ip and a — > ifi, assuring the formation of the d'Alembert operator 
when the timelikc and spacelike fluctuations are added with equal weight (this 

1 This transformation is similar to the continuation p —* ip' in nonrelativistic tunnelling, 
for which the analytic continuation appears as an instanton. 
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equality is consistent with the natural assumption, in this case, of an equal dis- 
tribution between spacelike and timelike contributions). The preservation of the 
magnitude of the interval reflects the conservation of a mass-like property which 
remains, as an intrinsic property of the particle, for both spacelike and timelike 
jumps. As mentioned before, one recalls the role of analytic continuation in 
quantum field theory; for the well known Wick rotation [e.g. H51], however, in 
that case, only the 0-component is analytically continued and no clear direct 
physical idea or quantity is associated with it. In the RBP the identification 
of the imaginary 4-momentum is dynamical in origin. It is due to the Lorentz 
structure of spacetime, which distinguishes the transitions > 1 from those 
with < 1. Though one may object to the association of Ax 1 * with a dy- 
namical momentum (since the instantaneous derivative is not defined for 
a Brownian process) the Brownian motion could be understood as an approx- 
imation to a microscopic process, just as it appears in Einstein's famous work 
in 1905 yQ, where it is assumed that the Brownian motion is produced by colli- 
sions. The effective conservation of E 2 — p 2 as a real quantity in both timelike 
and spacelike processes suggests that it is a physical property which preserves 
its meaning in both sectors. 

With these assumptions, the cross-term in hyperbolic functions cancels in 
the sum, which now takes the form 

Taking into account the drift term in Eq. one then finds the relativistic 
Fokker-Planck equation 

dp(x,r) , d . ,. d d , , . 

^-l-^F + W^e^**'^ (14) 

where djdx^ operates on both /3 M and p. 

We see that the procedure we have followed, identifying a 2 = —p 2 and as- 
suming equal weight, permits us to construct the Lorentz invariant 
d'Alcmbertian operator, as required for obtaining a relativistically covariant 
diffusion equation. 

To see this process in terms of a higher symmetry, let us define the invariant 
k 2 = Ef — p 2 > for the timelike case; our requirement is then that E 2 — p 2 = 
— k 2 for the spacelike case. In the framework of a larger group that includes 
k as part of a three vector (E,K,p), the relation for the timelike case can be 
considered in terms of the invariant of the subgroup Oil, 2), i.e., E 2 — k 2 — p 2 . 
The change in sign for the spacelike case yields the invariant E 2 + k 2 — p 2 ; we 
designate the corresponding symmetry (keeping the order of E andp) as 0(2, 1). 
These two groups may be thought of as subgroups of 0(2, 2) , where there exists a 
transformation which changes the sign of the metric of the subgroups holding the 
quantity k 2 constant. The kinematic constraints we have imposed correspond to 
setting these invariants to zero (the zero interval in the 2 + 1 and 1 + 2 spaces) . 
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The constraint we have placed on the relation of the timelike and spacclikc 
invariants derives from the properties of the distribution function and the re- 
quirement of obtaining the d'Alembert operator, i.e, Lorentz covariance of the 
diffusion equation. It appears that in order for the Brownian motion to result 
in a covariant diffusion equation, the distribution function has a higher sym- 
metry reflecting the necessary constraints. The transformations E — > iE' and 
p — > ip' used above would then correspond to analytic continuations from one 
(subgroup) sector to another. We shall see a similar structure in the 3 + 1 case, 
where the groups involved can be identified with the symmetries of the U(l) 
gauge fields associated with the quantum Stueckelberg-Schrodinger equation. 



5 The Invariance Properties of the Process 



Before formulating the 3+1 dimensional Brownian process, let us investigate the 
Lorentz invariance of the process and the correlation functions. The averaging 
operations are summations with weights (probability) assigned to each quantity 
in the sum. The sums in the continuum are, of course, expressed by integrals. 
If we wish to assign a relativistic Gaussian distribution function then the hy- 
perbolic angle integration is infinite unless we introduce a cutoff. The question 
then arises whether our process is invariant or not. 

We will show that we can use an arbitrary non-invariant (scalar) probabil- 
ity distribution (for example, a cutoff on the hyperbolic angle) and still obtain 
Lorentz invariant averages, using the imaginary representations of the 'unphysi- 
cal jumps'. For example, < Aw M > stands for a summation with a scalar weight 
(given by the density) over all the vectors Aw M , in the domain. It is therefore 
a vector. Moreover under the imaginary representation of spacelike increments 
relative to the timelike ones (here we assume the timelike jumps physical), Aw' 1 
is a simple vector function over all spacetime which has the following form: 

Aw' 1 = Aw' M , Aw'' 1 timelike 

Aw' 1 = iAw'' 1 , AW spacelike. (15) 

where the Aw /At are real. The quantity < Aw M > (formally written as a discrete 
sum) is given therefore by: 

< Awf > = J2 p ( Aw ) Aw,1 = P'(Aw')Aw ,f * +i P'{Aw')Aw 

timelike spacclikc 

= < Aw'^ >timclikc +i < Aw' M > sp acclikc 

where P(Aw) (or P'(Aw')) is the probability(weight) of having the vector Aw 
(or Aw'). The two vectors in the last equality in Eq. Ijl6|l are just normal 
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Lorentz vectors. If we now pick a distribution in a given frame for which the 
average of each of them (independent of the other) is zero then < Aw^ >= is 
true in all frames since the 0- vector is Lorentz invariant. 

Building the second correlation, with the assumption of no correlation be- 
tween spacelike jumps and timelike jumps, we find: 

< A^Aw" > = ^2P{Aw)Aw"Aw u = P\&w')Aw' tl Aw' u + 

time like 

+ i 2 P'{^w')Aw' tl Aw llJ 
spacelike 

— a timelike °spacelike ^ ^ (1?) 

where is the correlation tensor in each case. From the definition of Aiz/ M 
(a four vector) it follows that are real Lorentz tensors. The last equality in 
Eq. I|17|) is a demand that could be achieved for the general l+n case, by assum- 
ing that in a given frame there is an invariant Gaussian distribution where the 
distribution is uniform in all angles and that there is a cutoff in the hyperbolic 
angle. The sum of the two covariant tensors (each a result of summation on 
different sectors) is a Lorentz invariant tensor. The higher correlation functions 
do not interest us since they are of higher order in p and therefore in At and 
do not contribute to the Fokker-Plank equation. 



6 The Notion of 'Jumps' Versus a Continuous 
Process 

The mapping given in Eq. I|15|) leads necessarily to a deviation from the standard 
mathematical formulation of Brownian motion. There the probability that a 
particle starting at x at time r ending at x' at time r' is equal to the probability 
that the particle starts at x at time r passing through any possible intermediate 
point x" at time t" < r 1 and going from there to the point x' at time r' . This 
property is expressed in the Chapman-Kolmagorov equation [e.g. 

p(x, r, x', At') = f p(x, t, x", t") P (x", t" , x', t')c1 4 x", (18) 
Jn 

In the relativistic formulation the vector Aw' = x — x' could be a timelike 
vector therefore resulting in a real valued vector Aiu according to the mapping in 
Eq. (|15fl However, the two intermediate vectors Aw[ — x — x" and Aw' 2 = x" — x 
could be spacelike, and take the event out of the real manifold into a complex 
valued coordinate. In this case the Chapman-Kolmagorov equation does not 
hold, and the event may be found outside of the real manifold. In order to build a 
consistent process one must adopt the concept of 'Brownian jumps' which could 
be a result for example of a process in which the event (similar to Einstein's 
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original construction) undergoes collisions and for each collision, or 'jump' the 
mapping in Eq. (|15f) holds. Therefore at each point in the physical manifold 
the event may take any increment spacelike or timelike (with a possible complex 
valued contribution to the averages). However, although the vector leading from 
the initial point , say O, to the end point, A, may be spacelike and therefore 
be represented as an imaginary vector it is understood that the event arrives 
at the real spacetime point A, never leaving the real spacetime. This structure 
separates the two manifolds, spacetime which is real and represents the physical 
coordinates of the event and a complex space representing the processes the 
event undergoes going from one point to another. This structure differs in that 
sense from the mathematical formulation due to Wiener and others, but still it 
can be shown that the process is invariant on the average under decomposition 
into shorter time subprocesses . In other words, we consider the event starting 
at some arbitrary point, and going for some time At. We next decompose the 
time interval into M intervals, so that: 

M 

At, = At 

i=i 

We consider then the expression appearing in the Fokkcr-Plank equation 

< Ax^Ax" >: 

M M 

< Ax»Ax" >=< C£ A <) (E A <) > ( 19 ) 

»=i j=i 

Since we assume that any two non-equal time jumps are not correlated, i.e. 

< AxiAxj = > for i =/= j, which leaves only the equal time averages in the 
sum, 

M M 

< Ax^Ax" >= J2 < A < A < >= ^V" E Ar * = ct V"At (20) 

i=l i=l 

where a 2 is the diffusion constant and we used Eq. I|17|l going from the second 
piece in the equality to the third 

The notion of 'jumps' stimulates the consideration of discrete processes, 
which can also be formulated within the relativistic framework and leads, under 
certain assumptions, to a covariant Fokker-Plank equation. For example let us 
assume a physical process in which the 'jumps' occur in a very ordered way every 
tj seconds with a very small time spread (i.e. a very small probability that a 
collision occurs within a time different significantly from tj). Then, averaging 
the 'jumps' over a period r >> tj leads to: 

Nr., < t < (N + 1) TJ (21) 

This result is due to the fact that under our assumptions during the time r, N 
single 'jumps' within separation of each other of tj occurred. The average in 
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Eq. (|21l) does not change when r changes in less then tj; however if t,j is small 
then one can replace Eq. 121|) . with, 

< Am^Aw" >= o 2 rf v T 

(22) 

Therefore we recover the standard result for Brownian motion. However there 
is one very important difference which is the fact that r can be taken to be 
finitely small where in the standard Brownian process r can be actually taken 
to zero. This implies that higher order derivative terms enter into the resulting 
'diffusion' equation. For example for an isotropic homogeneous Gaussian dis- 
tribution there will be additional even order derivative operators beyond the 
second order (d'Alembert) with coefficients — 1) where n is the (even) 

order of the differential operator. Since both tj and a 2 are small these opera- 
tors could be neglected in general, though there might be special configurations 
in which their effect may be significant. In the following we assume that the 
tj is very small compared with the macroscopic scale and that the 'jumps' are 
practically ordered with zero spread, thus the approximation in Eq. ill' 21) is valid 
and no higher order terms are considered. 



7 Brownian motion in 3 + 1 dimensions 



In the 3 + 1 case, we again separate the jumps into timelike and spacelike types. 
The spacelike jumps may be parameterized, in a given frame, by 

Aw° — [i sinh a 

Aw 1 — /i cosh a cos <fr sin i? 

Aw 2 — [i cosh a sin <f> sin $ 

Aw 3 = /I cosh a cos 1} (23) 

We assume the four variables [A, a, 1?, <fi are independent random variables. 
In addition we demand in this frame that z9 and (j) uniformly distributed in 
their ranges (0,7r) and (0, 2tt), respectively. In this case, we may average over 
the trigonometric angles, i.e., "d and 4> and find that: 

2 

< Aw l2 >0, tf = < Aw 22 Au; 32 >0,^= ^-cosh 2 a 

o 

< Aid" 2 > <t> ^ = /i 2 sinh 2 a (24) 

We may obtain the averages over the trigonometric angles of the timelike jumps 
by replacing everywhere in Eq. 124|) 

cosh 2 a «-> sinh 2 a 
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to obtain 



< Au> l2 ><p^ = < Aw 22 >0 )1 j=< Aw 32 >(f>^= ^-sinh 2 a 

< Aw° 2 = cr 2 cosh 2 a, (25) 

where a is a real random variable, the invariant timelike interval. Assuming, as 
in the 1 + 1 case, that the likelihood of the jumps being in either the spacelike 
or (virtual) timelike phases are equal, and making an analytic continuation for 
which a 2 — > —A 2 , the total average of the operator O, including the contribu- 
tions of the remaining degrees of freedom /i, A and a is 

d 2 

< O > = (< u 2 >< sinh 2 a > - < A 2 >< cosh 2 a >) — ^ + 

-(< >< cosh 2 a > - < A 2 >< sinh 2 a >)A (26) 
o 

If we now insist that the operator < O > is invariant under Lorentz transfor- 
mations (i.e. the d'Alembertian) we impose the condition 

< ^ 2 >< sinh 2 a > - < A 2 >< cosh 2 a >= 
— o (< M 2 >< cosh 2 a > - < A 2 >< sinh 2 a >) (27) 

Using the fact that < cosh 2 a > — < sinh 2 a >= 1, and defining 
7 =< sinh 2 a >, we find that 

<A 2 >=i±i2<, 2 > (28) 

The Fokker-Planck equation then takes on the same form as in the 1 + 1 case, 
i.e., the form Eq. (|14J) . We remark that for the 1 + 1 case, one finds in the 
corresponding expression that the 3 in the denominator is replaced by unity, 
and the coefficients 4 are replaced by 2; in this case the requirement reduces to 
< /i 2 >=< A 2 > and there is no 7 dependence. 

We see that in the limit of a uniform distribution in a, for which 7 — > 00, 

< A 2 >^< /i 2 > . (29) 

In this case, the relativistic generalization of nonrclativistic Gaussian distribu- 

tion of the form e « is e ^, which is Lorentz invariant. 

The limiting case 7^0 corresponds to a stochastic process in which in 
the spacelike case there are no fluctuations in time, i.e., the process is that of 
a nonrelativistic Brownian motion. For the timelike case (recall that we have 
assumed the same distribution function over the hyperbolic variable) this limit 
implies that the fluctuations are entirely in the time direction. The limit 7 — > 00 
is Lorentz invariant, but the limit 7 — > can clearly be true only in a particular 
frame. 
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8 The Markov Relation and the 4D Gaussian 
Process 



In developing the previous ideas leading to the formulation of a RBP, we as- 
sumed that the probability distribution is consistent with the Markov property 
expressed in the Chapman-Kolmagorov equation [e.g. However, for the 

relativistic Gaussian it is not clear whether Eq. (|18|l holds. Therefore we now 
consider an alternative process, using the ideas developed above, resulting even- 
tually in the Klein-Gordon equation. Let us consider a 2D Gaussian process 
generated by a distribution of the form : 

P{W > dT) = 2^mr eX1P[ 2Ddr } (30) 

This distribution corresponds to a Markov process, a standard normalized Wiener 
process, where D is the diffusion constant. We now use the coordinate repre- 
sentation given in Eq. Q and Eq. JSJ for the timelike and spacelike sectors to 
transform the distribution function in Eq. (|30|) in both sectors to (use /i 2 in 
both cases): 

1 . — u 2 cosh 2a. . , 

:™p( ) (31) 



2-nDdr v 2Ddr 
where timelike 'jumps' are physical and the measure for both sectors is /id[ida 
Then, following Section 5, using Eq. (TSJ, we get for the combination of the 
timelike and spacelike contributions (with the appropriate sign) of the averages, 
say, A w 2 , 

. 2 1 f°° r 3 / -/i 2 cosh2a w , 
< Aw o > = TT-t^j- / / A* e M W da = 



where we integrated over /i first, using 



2nDdT J ^ 2Ddr 

— Ddr (32) 



H n exp{-a f i 2 ) = (33) 



and then integrated over a using 

J-co cosh 2 2a ^ ^ 

In a similar way one finds that (using Eq. (|15|) leading to the negative sign) 

< Awi 2 >= --Ddr (35) 

Since the probability distribution Eq. (|30|) is symmetric in Aw; in each sector 
< AwoAwi >— as well as the first moments. Therefore we get in this partic- 
ular frame a d'Alembertian. However, following the results of Section 5 we see 
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that it is an invariant result in all Lorentz frames (though in other frames the 
distribution may not appear to be Gaussian). 

Next we consider the application of the 4D form of Eq. (|3U[1 

1 -Aw 2 - Awi 2 - Aw 2 2 - Aw 3 2 
P(v> d -) = ^WW? eXP{ 2Dd? } (36) 

with measure /i 3 d^icosh 2 asin9d9dad<p for the spacelike sector and /i 3 d^isinh 2 asinOdQdadtfi 
for the timelike sector. 

However, now calculating < Awo 2 > for the timelike case, after averaging 
over the spatial angles 9 and <f> we find, using Eq. (|23[1 . 

. 2 1 f°° f°° 5 ,-M 2 cosh2a, l2 . l2 , , 

< Awo >= / / u exp( ; 1 cosh a sum adfida 

nD 2 (dT) 2 J J_ 0o 2Ddr 

(37) 

and for the spacelike case we get the same result since the spacelike parametriza- 
tion of Awo 2 is proportional to sinh 2 a and the spacelike volume element is pro- 
portional to cosh 2 a. Therefore if we use Eq. (|15|l . adding the contribution of 
the two sectors one obtains a complete cancellation to zero. In order to avoid 
this we extend Eq. (|15|l to the form 

Aw^ 1 = Aw'^ , Aw'' 1 timelike 

Aw^ 1 = i\Aw ,f * , Aw'' 1 spacelike (38) 

Before completing the calculation, we discuss the inclusion of the factor A, 
in Eq. I|38(l. Let us consider a classical (i.e. non-stochastic) event with a given 
value m 2 = Aw (i Ato m , moving in a timelike direction. It then changes its state 
of motion and starts moving in a spacelike direction; according to Eq. I|38|) 
m 2 changes into A 2 m 2 . Moreover, though the event may move according to a 
Gaussian distribution which makes no distinction between timelike and spacelike 
motions, the outcome of this motion as represented by the Aw, in Eq. (|15|) : 
Eq. (|38|l does distinguish the two phases of motion. We shall see that a specific 
value of A is required for the realization of the Fokker-Planck equation. 

The w manifold is complex and it is a function of the motion on the real 
manifold w' . Our macroscopic (physical) equations are written on the real plane 
of the w manifold. One can then visualize the flow of an event in spacetime 
similar to a motion of a particle in a cloud chamber. There as the particle 
moves the gas condenses, therefore the particle leaves a track. The track itself 
is not the particle but a result of the actual motion of the particle and its 
interaction with the gas in the cloud chamber. The track in the cloud chamber 
is analogous to the complex representation we use for the 'jumps'. 

We calculate first the expectation < Awq 2 > which is the total expectation, 
summed over the timelike and spacelike sectors. Averaging over the spherical 
angles 9, ip we get using Eq. l|3*7|l and Eq. J2HJl, 

. 2 i-A 2 ^,,^ f°° f°° 5 / -M 2 cosh2a N , 2 . , a , , 

< Aw > = D (dr) / / ^°exp{— ) cosh a smh aduda = 

TT Vn ./_™ 2Ddr 
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8(1 - A 2 ) „ , f°° cosh 2 a sinh 2 a , 

= -Ddr / = da (39) 

n J-oo cosh 2a 

where Eq. was used in the integration leading to the last equality in 
Eq. (|2jJ). Using 

2 1 

cosh a — -(cosh2a + l) 

sinh 2 a = i(cosh2a-l) (40) 



in Eq. 1|39|) and integrating over a we get 



< Aw 2 >= (1 X2) Ddrl (41) 

7T 2 



where we used 



da n 



h = 



h ^ I ~^T^ = T (42) 



x cosh 2a 2 
da 7r 
_oo cosh 3 2a 4 



We now calculate the expectation of < Awi 2 >. Averaging over the spherical 
angles 9, ip we get, using Eq. (|3T|) and Eq. (J3SJ, 

< Awi 2 > = — D 2 (dr) 2 f f [i 5 exp( - 2( * ) (sinh 4 a - A 2 cosh 4 a)dnda 
3?r J J_ oq 2Ddr 



(sinh 4 a — A 2 cosh 4 a) 
3vr J cosh 3 2a 



Ddr I 1 ~ 3n ' da (43) 



Using Eq. ijiOjl. Eq. (3^1. Eq. We get after integration over a, 

< Awi 2 >= ^i[(l - A 2 )^ - 4(1 + A 2 )] (44) 

In order to obtain the d'Alembertian we insist that < Awi 2 >= — < Awo 2 >, 
which leads to 

A 2 = *LZ± (45) 

3tt + 4 V 1 

Finally, substituting (for example) Eq. ijlSl in Eq. <|IT|) We find 

< Aw"Aw v >= if v AD dr = rf v Ddr (46) 

' 37T + 4 ' v ' 

where D, is the actual effective diffusion constant defined by 

t> S -*»- (47) 

3tt + 4 V 1 
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9 Discussion and Conclusions 



We have constructed a relativistic generalization of Brownian motion, using 
the invariant world-time, r, to order the Brownian fluctuations, and separated 
consideration of spacelike and timclike jumps to avoid the problems of negative 
second moments which might otherwise follow from the Minkowski signature. 
Associating the Brownian fluctuations with an underlying dynamical process, 
one may think of 7 discussed in the 3 + 1 case as an order parameter, where the 
distribution function (over a), associated with the velocities, is determined by 
the temperature of the underlying dynamical system (the result for the 1 + 1 case 
is independent of the distribution on the hyperbolic variable) . More generally it 
is suggestive to consider the possible thermodynamical effects of the 'medium' 
generating the relativistic Brownian fluctuations, following similar steps taken 
by Einstein ,1 ; in his famous work and verify whether any physical effect can be 
predicted. 

At equilibrium, where dp/dr — 0, the resulting diffusion equation turns into 
a classical wave equation which, in the absence of a drift term K^, is the wave 
equation for a massless field. An exponentially decreasing distribution in r of 
the form exp —kt would correspond to a Klein-Gordon equation for a particle 
in a tachyonic state (mass squared — k), for physical spacelike motion and for 
physical timclike motion to a particle with mass squared k. 

Choosing a cutoff in the hyperbolic angle, one finds a covariant moment and 
therefore, covariant differential operators. However the underlying process is 
not invariant, thus one can think of a special frame in which the hyperangular 
distribution is uniformly distributed around 0. Boosting breaks the symmetry 
of the hyperangular distribution, but since the the averages are tensor quantities 
the invariance properties are conserved, and therefore the Fokker-Plank equation 
(leading to the quantum equation) is invariant. This property is also used to 
construct the 4D Gaussian process. 

It was shown that a (Euclidian) Gaussian process with an appropriate (weighted) 
complex representation for the timclike and spacelike random motions can be 
used to achieve the covariant quantum equation, with the assurance that it is 
Markovian, (it is a relativistic generalization of the Wiener process). This leads 
to a 'cloud chamber-like' picture in which the event as it evolves leaves a track 
(carries a real or an imaginary phase), which is a representation of the actual 
motion, distinguishing the timelike and spacelike motion. 

In the classical Stueckelberg theory the timelike (forward or backward) prop- 
agation is associated with the standard particle or antiparticle interpretation, 
where spacelike propagation is needed whenever one discuss classical pair cre- 
ation or annihilation (with continuous passage from forward to backward motion 
in time). This suggests that the spacelike process may be associated with the 
annihilation and creation of pairs. Moreover, though the resulting macroscopic 
equation(i.e. on the level of the Fokker-Planck equation) is local and causal 
in the spacetime variables, the underlying microscopic process(i.e., on the level 
of the Brownian fluctuations) is not. It is however local and causal in r even 
at the microscopic level. This non-locality in t 'microscopically' may lead to a 
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mechanism providing correlations for the entangled state system. 

Nelson has shown that non-relativistic Brownian motion can be associated 
with a Schrodinger equation. Equipped with the procedures we presented here, 
constructing relativistic Brownian motion, Nelson's methods can be general- 
ized. One then can construct relativistic equations of Schrodinger (Schrodinger- 
Stueckelberg) type. The eigenvalue equations for these relativistic forms are also 
Klein-Gordon type equations. Moreover one can also generalize the case where 
the fluctuations are not correlated in different directions into the case where cor- 
relations exist, as discussed by Nelson for three dimensional Riemannian spaces. 
In this case the resulting equation will be a quantum equation in a curved Rie- 
mannian spacetime; as we have pointed out, the eikonal approximation to the 
solutions of such an equation contains the geodesic motion of classical general 
relativity. The medium supporting the Brownian motion may be identified with 
an "ether" for which the problem of local Lorentz symmetry is solved. This 
study opens up several tracks of possible research. Nelson, discussing the E.P.R 
0] system confronted the fact that such a system may be either described by a 
non-local Markov process or a local non-Markov process. The Markov process is 
simple to implement but Nelson was disturbed by the introduction of non-local 
interaction. However, the non-Markov process is very difficult to apply. The 
RPB developed here may bridge the two possibilities since an ordered (causal) 
Markov process in r may appear to be a non-Markovian (or possibly non-local 
and certainly non causal) process in t. For example for the Gaussian process, 
looking for the probability of finding the event changing its spatial position Ax 
after At has passed, one may integrate Eq. (|30J) over all r. This results how- 
ever, in Aa^+At^ wn i cn i s n °t integrable and therefore can not be normalized. 
This however is not surprising, since the probability of finding the event in Ax 
after At is not well defined (there may be several values of Aa; for a given At). 
For example the number of particles, as in Stueckelberg original construction 
jS] depends on the trajectory through which the point is reached. Defining an 
appropriate one particle probability resulting from the initial process occurring 
in t demands a restriction of the sample space before integrating over r i.e., 
using the conditional probability restricted to processes for which the event's t 
coordinate is monotonic in r (no pairs are created). 

Finally we would like to point out that generating a covariant quantum equa- 
tion through an RBP leads to a possible relation between quantum mechanics 
and gravitation. In the context of this work, the metric of gravity can appear 
as an anisotropy in the correlations that lead to quantum equations for which 
the ray, or eikonal approximation, corresponds to the classical geodesic flow of 
general relativity. It furthermore appears interesting to generalize Einstein's fa- 
mous work on this process introducing thermodynamic concepts to the resulting 
geometrical structure of the theory. 
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